\(\int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx\) [773]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (warning: unable to verify)
   Maple [C] (warning: unable to verify)
   Fricas [F]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 29, antiderivative size = 726 \[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=\frac {\sqrt {3+b} (c-d) \sqrt {c+d} (b c+15 d) E\left (\arcsin \left (\frac {\sqrt {3+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {3+b \sin (e+f x)}}\right )|\frac {(3-b) (c+d)}{(3+b) (c-d)}\right ) \sec (e+f x) \sqrt {-\frac {(b c-3 d) (1-\sin (e+f x))}{(c+d) (3+b \sin (e+f x))}} \sqrt {\frac {(b c-3 d) (1+\sin (e+f x))}{(c-d) (3+b \sin (e+f x))}} (3+b \sin (e+f x))}{4 (b c-3 d) d f}+\frac {\sqrt {c+d} \left (18 b c d+27 d^2-b^2 \left (c^2-4 d^2\right )\right ) \operatorname {EllipticPi}\left (\frac {b (c+d)}{(3+b) d},\arcsin \left (\frac {\sqrt {3+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {3+b \sin (e+f x)}}\right ),\frac {(3-b) (c+d)}{(3+b) (c-d)}\right ) \sec (e+f x) \sqrt {-\frac {(b c-3 d) (1-\sin (e+f x))}{(c+d) (3+b \sin (e+f x))}} \sqrt {\frac {(b c-3 d) (1+\sin (e+f x))}{(c-d) (3+b \sin (e+f x))}} (3+b \sin (e+f x))}{4 b \sqrt {3+b} d^2 f}-\frac {b (b c+15 d) \cos (e+f x) \sqrt {c+d \sin (e+f x)}}{4 d f \sqrt {3+b \sin (e+f x)}}-\frac {b \cos (e+f x) \sqrt {3+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}+\frac {(3+b)^{3/2} (9 d+b (c+2 d)) \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {3+b \sin (e+f x)}}{\sqrt {3+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(3+b) (c-d)}{(3-b) (c+d)}\right ) \sec (e+f x) \sqrt {\frac {(b c-3 d) (1-\sin (e+f x))}{(3+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-3 d) (1+\sin (e+f x))}{(3-b) (c+d \sin (e+f x))}} (c+d \sin (e+f x))}{4 b d \sqrt {c+d} f} \]

[Out]

1/4*(6*a*b*c*d+3*a^2*d^2-b^2*(c^2-4*d^2))*EllipticPi((a+b)^(1/2)*(c+d*sin(f*x+e))^(1/2)/(c+d)^(1/2)/(a+b*sin(f
*x+e))^(1/2),b*(c+d)/(a+b)/d,((a-b)*(c+d)/(a+b)/(c-d))^(1/2))*sec(f*x+e)*(a+b*sin(f*x+e))*(c+d)^(1/2)*(-(-a*d+
b*c)*(1-sin(f*x+e))/(c+d)/(a+b*sin(f*x+e)))^(1/2)*((-a*d+b*c)*(1+sin(f*x+e))/(c-d)/(a+b*sin(f*x+e)))^(1/2)/b/d
^2/f/(a+b)^(1/2)+1/4*(c-d)*(5*a*d+b*c)*EllipticE((a+b)^(1/2)*(c+d*sin(f*x+e))^(1/2)/(c+d)^(1/2)/(a+b*sin(f*x+e
))^(1/2),((a-b)*(c+d)/(a+b)/(c-d))^(1/2))*sec(f*x+e)*(a+b*sin(f*x+e))*(a+b)^(1/2)*(c+d)^(1/2)*(-(-a*d+b*c)*(1-
sin(f*x+e))/(c+d)/(a+b*sin(f*x+e)))^(1/2)*((-a*d+b*c)*(1+sin(f*x+e))/(c-d)/(a+b*sin(f*x+e)))^(1/2)/d/(-a*d+b*c
)/f+1/4*(a+b)^(3/2)*(3*a*d+b*(c+2*d))*EllipticF((c+d)^(1/2)*(a+b*sin(f*x+e))^(1/2)/(a+b)^(1/2)/(c+d*sin(f*x+e)
)^(1/2),((a+b)*(c-d)/(a-b)/(c+d))^(1/2))*sec(f*x+e)*(c+d*sin(f*x+e))*((-a*d+b*c)*(1-sin(f*x+e))/(a+b)/(c+d*sin
(f*x+e)))^(1/2)*(-(-a*d+b*c)*(1+sin(f*x+e))/(a-b)/(c+d*sin(f*x+e)))^(1/2)/b/d/f/(c+d)^(1/2)-1/4*b*(5*a*d+b*c)*
cos(f*x+e)*(c+d*sin(f*x+e))^(1/2)/d/f/(a+b*sin(f*x+e))^(1/2)-1/2*b*cos(f*x+e)*(a+b*sin(f*x+e))^(1/2)*(c+d*sin(
f*x+e))^(1/2)/f

Rubi [A] (verified)

Time = 1.62 (sec) , antiderivative size = 740, normalized size of antiderivative = 1.02, number of steps used = 7, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.241, Rules used = {2900, 3140, 3132, 2890, 3077, 2897, 3075} \[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=\frac {\sqrt {c+d} \left (3 a^2 d^2+6 a b c d-\left (b^2 \left (c^2-4 d^2\right )\right )\right ) \sec (e+f x) (a+b \sin (e+f x)) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (\sin (e+f x)+1)}{(c-d) (a+b \sin (e+f x))}} \operatorname {EllipticPi}\left (\frac {b (c+d)}{(a+b) d},\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right ),\frac {(a-b) (c+d)}{(a+b) (c-d)}\right )}{4 b d^2 f \sqrt {a+b}}+\frac {(a+b)^{3/2} (3 a d+b (c+2 d)) \sec (e+f x) (c+d \sin (e+f x)) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (\sin (e+f x)+1)}{(a-b) (c+d \sin (e+f x))}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right )}{4 b d f \sqrt {c+d}}+\frac {\sqrt {a+b} (c-d) \sqrt {c+d} (5 a d+b c) \sec (e+f x) (a+b \sin (e+f x)) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (\sin (e+f x)+1)}{(c-d) (a+b \sin (e+f x))}} E\left (\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right )|\frac {(a-b) (c+d)}{(a+b) (c-d)}\right )}{4 d f (b c-a d)}-\frac {b \cos (e+f x) \sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}-\frac {b (5 a d+b c) \cos (e+f x) \sqrt {c+d \sin (e+f x)}}{4 d f \sqrt {a+b \sin (e+f x)}} \]

[In]

Int[(a + b*Sin[e + f*x])^(3/2)*Sqrt[c + d*Sin[e + f*x]],x]

[Out]

(Sqrt[a + b]*(c - d)*Sqrt[c + d]*(b*c + 5*a*d)*EllipticE[ArcSin[(Sqrt[a + b]*Sqrt[c + d*Sin[e + f*x]])/(Sqrt[c
 + d]*Sqrt[a + b*Sin[e + f*x]])], ((a - b)*(c + d))/((a + b)*(c - d))]*Sec[e + f*x]*Sqrt[-(((b*c - a*d)*(1 - S
in[e + f*x]))/((c + d)*(a + b*Sin[e + f*x])))]*Sqrt[((b*c - a*d)*(1 + Sin[e + f*x]))/((c - d)*(a + b*Sin[e + f
*x]))]*(a + b*Sin[e + f*x]))/(4*d*(b*c - a*d)*f) + (Sqrt[c + d]*(6*a*b*c*d + 3*a^2*d^2 - b^2*(c^2 - 4*d^2))*El
lipticPi[(b*(c + d))/((a + b)*d), ArcSin[(Sqrt[a + b]*Sqrt[c + d*Sin[e + f*x]])/(Sqrt[c + d]*Sqrt[a + b*Sin[e
+ f*x]])], ((a - b)*(c + d))/((a + b)*(c - d))]*Sec[e + f*x]*Sqrt[-(((b*c - a*d)*(1 - Sin[e + f*x]))/((c + d)*
(a + b*Sin[e + f*x])))]*Sqrt[((b*c - a*d)*(1 + Sin[e + f*x]))/((c - d)*(a + b*Sin[e + f*x]))]*(a + b*Sin[e + f
*x]))/(4*b*Sqrt[a + b]*d^2*f) - (b*(b*c + 5*a*d)*Cos[e + f*x]*Sqrt[c + d*Sin[e + f*x]])/(4*d*f*Sqrt[a + b*Sin[
e + f*x]]) - (b*Cos[e + f*x]*Sqrt[a + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]])/(2*f) + ((a + b)^(3/2)*(3*a*d
+ b*(c + 2*d))*EllipticF[ArcSin[(Sqrt[c + d]*Sqrt[a + b*Sin[e + f*x]])/(Sqrt[a + b]*Sqrt[c + d*Sin[e + f*x]])]
, ((a + b)*(c - d))/((a - b)*(c + d))]*Sec[e + f*x]*Sqrt[((b*c - a*d)*(1 - Sin[e + f*x]))/((a + b)*(c + d*Sin[
e + f*x]))]*Sqrt[-(((b*c - a*d)*(1 + Sin[e + f*x]))/((a - b)*(c + d*Sin[e + f*x])))]*(c + d*Sin[e + f*x]))/(4*
b*d*Sqrt[c + d]*f)

Rule 2890

Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]/Sqrt[(c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]], x_Symbol] :> Simp[
2*((a + b*Sin[e + f*x])/(d*f*Rt[(a + b)/(c + d), 2]*Cos[e + f*x]))*Sqrt[(b*c - a*d)*((1 + Sin[e + f*x])/((c -
d)*(a + b*Sin[e + f*x])))]*Sqrt[(-(b*c - a*d))*((1 - Sin[e + f*x])/((c + d)*(a + b*Sin[e + f*x])))]*EllipticPi
[b*((c + d)/(d*(a + b))), ArcSin[Rt[(a + b)/(c + d), 2]*(Sqrt[c + d*Sin[e + f*x]]/Sqrt[a + b*Sin[e + f*x]])],
(a - b)*((c + d)/((a + b)*(c - d)))], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2
, 0] && NeQ[c^2 - d^2, 0] && PosQ[(a + b)/(c + d)]

Rule 2897

Int[1/(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Si
mp[2*((c + d*Sin[e + f*x])/(f*(b*c - a*d)*Rt[(c + d)/(a + b), 2]*Cos[e + f*x]))*Sqrt[(b*c - a*d)*((1 - Sin[e +
 f*x])/((a + b)*(c + d*Sin[e + f*x])))]*Sqrt[(-(b*c - a*d))*((1 + Sin[e + f*x])/((a - b)*(c + d*Sin[e + f*x]))
)]*EllipticF[ArcSin[Rt[(c + d)/(a + b), 2]*(Sqrt[a + b*Sin[e + f*x]]/Sqrt[c + d*Sin[e + f*x]])], (a + b)*((c -
 d)/((a - b)*(c + d)))], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0] && NeQ[c
^2 - d^2, 0] && PosQ[(c + d)/(a + b)]

Rule 2900

Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Si
mp[(-b)*Cos[e + f*x]*(a + b*Sin[e + f*x])^(m - 1)*((c + d*Sin[e + f*x])^n/(f*(m + n))), x] + Dist[1/(d*(m + n)
), Int[(a + b*Sin[e + f*x])^(m - 2)*(c + d*Sin[e + f*x])^(n - 1)*Simp[a^2*c*d*(m + n) + b*d*(b*c*(m - 1) + a*d
*n) + (a*d*(2*b*c + a*d)*(m + n) - b*d*(a*c - b*d*(m + n - 1)))*Sin[e + f*x] + b*d*(b*c*n + a*d*(2*m + n - 1))
*Sin[e + f*x]^2, x], x], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0] && NeQ[c
^2 - d^2, 0] && LtQ[0, m, 2] && LtQ[-1, n, 2] && NeQ[m + n, 0] && (IntegerQ[m] || IntegersQ[2*m, 2*n])

Rule 3075

Int[((A_) + (B_.)*sin[(e_.) + (f_.)*(x_)])/(((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(3/2)*Sqrt[(c_) + (d_.)*sin
[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Simp[-2*A*(c - d)*((a + b*Sin[e + f*x])/(f*(b*c - a*d)^2*Rt[(a + b)/(c +
d), 2]*Cos[e + f*x]))*Sqrt[(b*c - a*d)*((1 + Sin[e + f*x])/((c - d)*(a + b*Sin[e + f*x])))]*Sqrt[(-(b*c - a*d)
)*((1 - Sin[e + f*x])/((c + d)*(a + b*Sin[e + f*x])))]*EllipticE[ArcSin[Rt[(a + b)/(c + d), 2]*(Sqrt[c + d*Sin
[e + f*x]]/Sqrt[a + b*Sin[e + f*x]])], (a - b)*((c + d)/((a + b)*(c - d)))], x] /; FreeQ[{a, b, c, d, e, f, A,
 B}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && EqQ[A, B] && PosQ[(a + b)/(c + d)]

Rule 3077

Int[((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])/(((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(3/2)*Sqrt[(c_) + (d_.)*s
in[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Dist[(A - B)/(a - b), Int[1/(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e
+ f*x]]), x], x] - Dist[(A*b - a*B)/(a - b), Int[(1 + Sin[e + f*x])/((a + b*Sin[e + f*x])^(3/2)*Sqrt[c + d*Sin
[e + f*x]]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0] && NeQ[c^2
 - d^2, 0] && NeQ[A, B]

Rule 3132

Int[((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) + (f_.)*(x_)]^2)/(((a_.) + (b_.)*sin[(e_.) + (f_.
)*(x_)])^(3/2)*Sqrt[(c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Dist[C/b^2, Int[Sqrt[a + b*Sin[e + f
*x]]/Sqrt[c + d*Sin[e + f*x]], x], x] + Dist[1/b^2, Int[(A*b^2 - a^2*C + b*(b*B - 2*a*C)*Sin[e + f*x])/((a + b
*Sin[e + f*x])^(3/2)*Sqrt[c + d*Sin[e + f*x]]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C}, x] && NeQ[b*c - a
*d, 0] && NeQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0]

Rule 3140

Int[((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) + (f_.)*(x_)]^2)/(Sqrt[(a_.) + (b_.)*sin[(e_.) +
(f_.)*(x_)]]*Sqrt[(c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Simp[(-C)*Cos[e + f*x]*(Sqrt[c + d*Sin[
e + f*x]]/(d*f*Sqrt[a + b*Sin[e + f*x]])), x] + Dist[1/(2*d), Int[(1/((a + b*Sin[e + f*x])^(3/2)*Sqrt[c + d*Si
n[e + f*x]]))*Simp[2*a*A*d - C*(b*c - a*d) - 2*(a*c*C - d*(A*b + a*B))*Sin[e + f*x] + (2*b*B*d - C*(b*c + a*d)
)*Sin[e + f*x]^2, x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0
] && NeQ[c^2 - d^2, 0]

Rubi steps \begin{align*} \text {integral}& = -\frac {b \cos (e+f x) \sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}+\frac {\int \frac {\frac {1}{2} d \left (4 a^2 c+b^2 c+a b d\right )+d \left (3 a b c+2 a^2 d+b^2 d\right ) \sin (e+f x)+\frac {1}{2} b d (b c+5 a d) \sin ^2(e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx}{2 d} \\ & = -\frac {b (b c+5 a d) \cos (e+f x) \sqrt {c+d \sin (e+f x)}}{4 d f \sqrt {a+b \sin (e+f x)}}-\frac {b \cos (e+f x) \sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}+\frac {\int \frac {-\frac {1}{2} d \left (b^3 c^2-8 a^3 c d+2 a b^2 c d-7 a^2 b d^2\right )+d \left (5 a^2 b c d+b^3 c d+4 a^3 d^2-a b^2 \left (c^2-3 d^2\right )\right ) \sin (e+f x)+\frac {1}{2} b d \left (6 a b c d+3 a^2 d^2-b^2 \left (c^2-4 d^2\right )\right ) \sin ^2(e+f x)}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}} \, dx}{4 d^2} \\ & = -\frac {b (b c+5 a d) \cos (e+f x) \sqrt {c+d \sin (e+f x)}}{4 d f \sqrt {a+b \sin (e+f x)}}-\frac {b \cos (e+f x) \sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}+\frac {\int \frac {-\frac {1}{2} b^2 d \left (b^3 c^2-8 a^3 c d+2 a b^2 c d-7 a^2 b d^2\right )-\frac {1}{2} a^2 b d \left (6 a b c d+3 a^2 d^2-b^2 \left (c^2-4 d^2\right )\right )+b \left (-a b d \left (6 a b c d+3 a^2 d^2-b^2 \left (c^2-4 d^2\right )\right )+b d \left (5 a^2 b c d+b^3 c d+4 a^3 d^2-a b^2 \left (c^2-3 d^2\right )\right )\right ) \sin (e+f x)}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}} \, dx}{4 b^2 d^2}+\frac {1}{8} \left (6 a c-\frac {b c^2}{d}+\frac {3 a^2 d}{b}+4 b d\right ) \int \frac {\sqrt {a+b \sin (e+f x)}}{\sqrt {c+d \sin (e+f x)}} \, dx \\ & = \frac {\sqrt {c+d} \left (6 a c-\frac {b c^2}{d}+\frac {3 a^2 d}{b}+4 b d\right ) \operatorname {EllipticPi}\left (\frac {b (c+d)}{(a+b) d},\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right ),\frac {(a-b) (c+d)}{(a+b) (c-d)}\right ) \sec (e+f x) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (1+\sin (e+f x))}{(c-d) (a+b \sin (e+f x))}} (a+b \sin (e+f x))}{4 \sqrt {a+b} d f}-\frac {b (b c+5 a d) \cos (e+f x) \sqrt {c+d \sin (e+f x)}}{4 d f \sqrt {a+b \sin (e+f x)}}-\frac {b \cos (e+f x) \sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}-\frac {((a+b) (b c-a d) (b c+5 a d)) \int \frac {1+\sin (e+f x)}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}} \, dx}{8 d}--\frac {\left (-\frac {1}{2} b^2 d \left (b^3 c^2-8 a^3 c d+2 a b^2 c d-7 a^2 b d^2\right )-\frac {1}{2} a^2 b d \left (6 a b c d+3 a^2 d^2-b^2 \left (c^2-4 d^2\right )\right )-b \left (-a b d \left (6 a b c d+3 a^2 d^2-b^2 \left (c^2-4 d^2\right )\right )+b d \left (5 a^2 b c d+b^3 c d+4 a^3 d^2-a b^2 \left (c^2-3 d^2\right )\right )\right )\right ) \int \frac {1}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx}{4 (a-b) b^2 d^2} \\ & = \frac {\sqrt {a+b} (c-d) \sqrt {c+d} (b c+5 a d) E\left (\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right )|\frac {(a-b) (c+d)}{(a+b) (c-d)}\right ) \sec (e+f x) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (1+\sin (e+f x))}{(c-d) (a+b \sin (e+f x))}} (a+b \sin (e+f x))}{4 d (b c-a d) f}+\frac {\sqrt {c+d} \left (6 a c-\frac {b c^2}{d}+\frac {3 a^2 d}{b}+4 b d\right ) \operatorname {EllipticPi}\left (\frac {b (c+d)}{(a+b) d},\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right ),\frac {(a-b) (c+d)}{(a+b) (c-d)}\right ) \sec (e+f x) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (1+\sin (e+f x))}{(c-d) (a+b \sin (e+f x))}} (a+b \sin (e+f x))}{4 \sqrt {a+b} d f}-\frac {b (b c+5 a d) \cos (e+f x) \sqrt {c+d \sin (e+f x)}}{4 d f \sqrt {a+b \sin (e+f x)}}-\frac {b \cos (e+f x) \sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}+\frac {(a+b)^{3/2} (3 a d+b (c+2 d)) \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right ) \sec (e+f x) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (1+\sin (e+f x))}{(a-b) (c+d \sin (e+f x))}} (c+d \sin (e+f x))}{4 b d \sqrt {c+d} f} \\ \end{align*}

Mathematica [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(1849\) vs. \(2(726)=1452\).

Time = 13.73 (sec) , antiderivative size = 1849, normalized size of antiderivative = 2.55 \[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=-\frac {b \cos (e+f x) \sqrt {3+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}{2 f}+\frac {\frac {4 (-b c+3 d) \left (-72 c-3 b^2 c-21 b d\right ) \sqrt {\frac {(c+d) \cot ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right )}{-c+d}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{\sqrt {2}}\right ),\frac {2 (-b c+3 d)}{(3+b) (-c+d)}\right ) \sec (e+f x) \sin ^4\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) \sqrt {\frac {(c+d) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (3+b \sin (e+f x))}{-b c+3 d}} \sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{(3+b) (c+d) \sqrt {3+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}+4 (-b c+3 d) \left (-36 b c-72 d-4 b^2 d\right ) \left (\frac {\sqrt {\frac {(c+d) \cot ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right )}{-c+d}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{\sqrt {2}}\right ),\frac {2 (-b c+3 d)}{(3+b) (-c+d)}\right ) \sec (e+f x) \sin ^4\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) \sqrt {\frac {(c+d) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (3+b \sin (e+f x))}{-b c+3 d}} \sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{(3+b) (c+d) \sqrt {3+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}-\frac {\sqrt {\frac {(c+d) \cot ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right )}{-c+d}} \operatorname {EllipticPi}\left (\frac {-b c+3 d}{(3+b) d},\arcsin \left (\frac {\sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{\sqrt {2}}\right ),\frac {2 (-b c+3 d)}{(3+b) (-c+d)}\right ) \sec (e+f x) \sin ^4\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) \sqrt {\frac {(c+d) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (3+b \sin (e+f x))}{-b c+3 d}} \sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{(3+b) d \sqrt {3+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )-2 \left (b^2 c+15 b d\right ) \left (\frac {\cos (e+f x) \sqrt {c+d \sin (e+f x)}}{d \sqrt {3+b \sin (e+f x)}}+\frac {\sqrt {\frac {3-b}{3+b}} (3+b) \cos \left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) E\left (\arcsin \left (\frac {\sqrt {\frac {3-b}{3+b}} \sin \left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right )}{\sqrt {\frac {3+b \sin (e+f x)}{3+b}}}\right )|\frac {2 (-b c+3 d)}{(3-b) (c+d)}\right ) \sqrt {c+d \sin (e+f x)}}{b d \sqrt {\frac {(3+b) \cos ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right )}{3+b \sin (e+f x)}} \sqrt {3+b \sin (e+f x)} \sqrt {\frac {3+b \sin (e+f x)}{3+b}} \sqrt {\frac {(3+b) (c+d \sin (e+f x))}{(c+d) (3+b \sin (e+f x))}}}-\frac {2 (-b c+3 d) \left (\frac {((3+b) c+3 d) \sqrt {\frac {(c+d) \cot ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right )}{-c+d}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{\sqrt {2}}\right ),\frac {2 (-b c+3 d)}{(3+b) (-c+d)}\right ) \sec (e+f x) \sin ^4\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) \sqrt {\frac {(c+d) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (3+b \sin (e+f x))}{-b c+3 d}} \sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{(3+b) (c+d) \sqrt {3+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}-\frac {(b c+3 d) \sqrt {\frac {(c+d) \cot ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right )}{-c+d}} \operatorname {EllipticPi}\left (\frac {-b c+3 d}{(3+b) d},\arcsin \left (\frac {\sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{\sqrt {2}}\right ),\frac {2 (-b c+3 d)}{(3+b) (-c+d)}\right ) \sec (e+f x) \sin ^4\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) \sqrt {\frac {(c+d) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (3+b \sin (e+f x))}{-b c+3 d}} \sqrt {\frac {(-3-b) \csc ^2\left (\frac {1}{2} \left (-e+\frac {\pi }{2}-f x\right )\right ) (c+d \sin (e+f x))}{-b c+3 d}}}{(3+b) d \sqrt {3+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{b d}\right )}{8 f} \]

[In]

Integrate[(3 + b*Sin[e + f*x])^(3/2)*Sqrt[c + d*Sin[e + f*x]],x]

[Out]

-1/2*(b*Cos[e + f*x]*Sqrt[3 + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]])/f + ((4*(-(b*c) + 3*d)*(-72*c - 3*b^2*
c - 21*b*d)*Sqrt[((c + d)*Cot[(-e + Pi/2 - f*x)/2]^2)/(-c + d)]*EllipticF[ArcSin[Sqrt[((-3 - b)*Csc[(-e + Pi/2
 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)]/Sqrt[2]], (2*(-(b*c) + 3*d))/((3 + b)*(-c + d))]*Sec[e + f*
x]*Sin[(-e + Pi/2 - f*x)/2]^4*Sqrt[((c + d)*Csc[(-e + Pi/2 - f*x)/2]^2*(3 + b*Sin[e + f*x]))/(-(b*c) + 3*d)]*S
qrt[((-3 - b)*Csc[(-e + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)])/((3 + b)*(c + d)*Sqrt[3 + b*Si
n[e + f*x]]*Sqrt[c + d*Sin[e + f*x]]) + 4*(-(b*c) + 3*d)*(-36*b*c - 72*d - 4*b^2*d)*((Sqrt[((c + d)*Cot[(-e +
Pi/2 - f*x)/2]^2)/(-c + d)]*EllipticF[ArcSin[Sqrt[((-3 - b)*Csc[(-e + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(
-(b*c) + 3*d)]/Sqrt[2]], (2*(-(b*c) + 3*d))/((3 + b)*(-c + d))]*Sec[e + f*x]*Sin[(-e + Pi/2 - f*x)/2]^4*Sqrt[(
(c + d)*Csc[(-e + Pi/2 - f*x)/2]^2*(3 + b*Sin[e + f*x]))/(-(b*c) + 3*d)]*Sqrt[((-3 - b)*Csc[(-e + Pi/2 - f*x)/
2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)])/((3 + b)*(c + d)*Sqrt[3 + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]]
) - (Sqrt[((c + d)*Cot[(-e + Pi/2 - f*x)/2]^2)/(-c + d)]*EllipticPi[(-(b*c) + 3*d)/((3 + b)*d), ArcSin[Sqrt[((
-3 - b)*Csc[(-e + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)]/Sqrt[2]], (2*(-(b*c) + 3*d))/((3 + b)
*(-c + d))]*Sec[e + f*x]*Sin[(-e + Pi/2 - f*x)/2]^4*Sqrt[((c + d)*Csc[(-e + Pi/2 - f*x)/2]^2*(3 + b*Sin[e + f*
x]))/(-(b*c) + 3*d)]*Sqrt[((-3 - b)*Csc[(-e + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)])/((3 + b)
*d*Sqrt[3 + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]])) - 2*(b^2*c + 15*b*d)*((Cos[e + f*x]*Sqrt[c + d*Sin[e +
f*x]])/(d*Sqrt[3 + b*Sin[e + f*x]]) + (Sqrt[(3 - b)/(3 + b)]*(3 + b)*Cos[(-e + Pi/2 - f*x)/2]*EllipticE[ArcSin
[(Sqrt[(3 - b)/(3 + b)]*Sin[(-e + Pi/2 - f*x)/2])/Sqrt[(3 + b*Sin[e + f*x])/(3 + b)]], (2*(-(b*c) + 3*d))/((3
- b)*(c + d))]*Sqrt[c + d*Sin[e + f*x]])/(b*d*Sqrt[((3 + b)*Cos[(-e + Pi/2 - f*x)/2]^2)/(3 + b*Sin[e + f*x])]*
Sqrt[3 + b*Sin[e + f*x]]*Sqrt[(3 + b*Sin[e + f*x])/(3 + b)]*Sqrt[((3 + b)*(c + d*Sin[e + f*x]))/((c + d)*(3 +
b*Sin[e + f*x]))]) - (2*(-(b*c) + 3*d)*((((3 + b)*c + 3*d)*Sqrt[((c + d)*Cot[(-e + Pi/2 - f*x)/2]^2)/(-c + d)]
*EllipticF[ArcSin[Sqrt[((-3 - b)*Csc[(-e + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)]/Sqrt[2]], (2
*(-(b*c) + 3*d))/((3 + b)*(-c + d))]*Sec[e + f*x]*Sin[(-e + Pi/2 - f*x)/2]^4*Sqrt[((c + d)*Csc[(-e + Pi/2 - f*
x)/2]^2*(3 + b*Sin[e + f*x]))/(-(b*c) + 3*d)]*Sqrt[((-3 - b)*Csc[(-e + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/
(-(b*c) + 3*d)])/((3 + b)*(c + d)*Sqrt[3 + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]]) - ((b*c + 3*d)*Sqrt[((c +
 d)*Cot[(-e + Pi/2 - f*x)/2]^2)/(-c + d)]*EllipticPi[(-(b*c) + 3*d)/((3 + b)*d), ArcSin[Sqrt[((-3 - b)*Csc[(-e
 + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)]/Sqrt[2]], (2*(-(b*c) + 3*d))/((3 + b)*(-c + d))]*Sec
[e + f*x]*Sin[(-e + Pi/2 - f*x)/2]^4*Sqrt[((c + d)*Csc[(-e + Pi/2 - f*x)/2]^2*(3 + b*Sin[e + f*x]))/(-(b*c) +
3*d)]*Sqrt[((-3 - b)*Csc[(-e + Pi/2 - f*x)/2]^2*(c + d*Sin[e + f*x]))/(-(b*c) + 3*d)])/((3 + b)*d*Sqrt[3 + b*S
in[e + f*x]]*Sqrt[c + d*Sin[e + f*x]])))/(b*d)))/(8*f)

Maple [C] (warning: unable to verify)

Result contains complex when optimal does not.

Time = 12.45 (sec) , antiderivative size = 250087, normalized size of antiderivative = 344.47

method result size
default \(\text {Expression too large to display}\) \(250087\)

[In]

int((a+b*sin(f*x+e))^(3/2)*(c+d*sin(f*x+e))^(1/2),x,method=_RETURNVERBOSE)

[Out]

result too large to display

Fricas [F]

\[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=\int { {\left (b \sin \left (f x + e\right ) + a\right )}^{\frac {3}{2}} \sqrt {d \sin \left (f x + e\right ) + c} \,d x } \]

[In]

integrate((a+b*sin(f*x+e))^(3/2)*(c+d*sin(f*x+e))^(1/2),x, algorithm="fricas")

[Out]

integral((b*sin(f*x + e) + a)^(3/2)*sqrt(d*sin(f*x + e) + c), x)

Sympy [F]

\[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=\int \left (a + b \sin {\left (e + f x \right )}\right )^{\frac {3}{2}} \sqrt {c + d \sin {\left (e + f x \right )}}\, dx \]

[In]

integrate((a+b*sin(f*x+e))**(3/2)*(c+d*sin(f*x+e))**(1/2),x)

[Out]

Integral((a + b*sin(e + f*x))**(3/2)*sqrt(c + d*sin(e + f*x)), x)

Maxima [F]

\[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=\int { {\left (b \sin \left (f x + e\right ) + a\right )}^{\frac {3}{2}} \sqrt {d \sin \left (f x + e\right ) + c} \,d x } \]

[In]

integrate((a+b*sin(f*x+e))^(3/2)*(c+d*sin(f*x+e))^(1/2),x, algorithm="maxima")

[Out]

integrate((b*sin(f*x + e) + a)^(3/2)*sqrt(d*sin(f*x + e) + c), x)

Giac [F]

\[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=\int { {\left (b \sin \left (f x + e\right ) + a\right )}^{\frac {3}{2}} \sqrt {d \sin \left (f x + e\right ) + c} \,d x } \]

[In]

integrate((a+b*sin(f*x+e))^(3/2)*(c+d*sin(f*x+e))^(1/2),x, algorithm="giac")

[Out]

integrate((b*sin(f*x + e) + a)^(3/2)*sqrt(d*sin(f*x + e) + c), x)

Mupad [F(-1)]

Timed out. \[ \int (3+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)} \, dx=\int {\left (a+b\,\sin \left (e+f\,x\right )\right )}^{3/2}\,\sqrt {c+d\,\sin \left (e+f\,x\right )} \,d x \]

[In]

int((a + b*sin(e + f*x))^(3/2)*(c + d*sin(e + f*x))^(1/2),x)

[Out]

int((a + b*sin(e + f*x))^(3/2)*(c + d*sin(e + f*x))^(1/2), x)